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aiAPTER I 
INTRODUCTION 

1.1 Motivation and Relevance of Repo-nt In April of 1974 the National 
Aeronautics and Space Administration conducted a high altitude balloon 
experiment called LACATE (Lower Atmosphere Composition and Tenperaiture 
Experimt'nt) which enployed an infrared radiometer to sense remotely 
vertical profiles of the concentrations of selected atmospheric trace 
constituents and tenperature. The constituents were measured by invert- 
ing infrared radiance profile of the earth's horieon. The radionveter line 
of sight was scanned vertically across the horizon at approximately 0,25° 
per second, requiring 30 seconds to acquire a coirplete radiance profile. 

The specifications required that the relative vertical position of 

the data points making up a profile be known to approximately 30 arc 
seconds. The general description of the balloon system for acoonplishing 
the mission is given in reference (1), refer figures 1.1-1 and 1.1-2. 

In order to fix the orientation of the line of sight of the 
radic»neter, it is necessary to be able to determine the configuration of 
the platform in space, i.e, the attitude of the system. This can be accom- 
plished by simulating the balloon system and using the gyro output in 
conjunction with a parameter estimation process. Simulation of the balloon 
system requires a nathematical model plus axulysis of the model. Ihe 
required mathematical model has already been developed for use with the 
system simulation process and the details are described in reference (1) 
pages 9-38. The attitude of the balloon system can be determined once tiie 
initial conditions; i.e, the initial state, is known, 

1.2 Objective of the Beport The main objective of this report will be to 
develop a process to determine the unknown initial state parameters by 
enploying the outpxit of the system mathematical model in conjunction with 
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the output obtained from the instrumentation system. This system consisted 
of three orthogonally oriented rate gyros and magnetometer which were 
fixed to the balloon platform. The development of this process involves 
2 major steps. First, a parameter determination process must be developed 
which solves for the un3<nown initial state parameters of the system which, 
in turn, give the best fit to the data obtained from the instrumentation 
system. Second, a simulation process must be developed for computing 
the state of the balloon system which can then be compared with the 
actual state of the system. 

The body of this report consists of 2 major parts. The first part 
discusses the attitude determination process in general. The details of 
this part include 

(i) System Simulation Discussion of the method of solution and 
computation of the output of interest (i.e. angular velocity ocmponents along 
the platform axis) . 

(xi) Canputation of System Natural Frequencies . Discussion of a 

method for solving system eigenvalues. 

(iii) Optimization Technique . Discussion of the Hooke and Jeeves 

direct seardi method. 

(iv) Verification Process . Discussion of a method for verification of 
the attitude determination process. 

In the second part of the report all the numerical data and results are 
presented and discussed. 
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CHAPn^R II 

ATTITUDE DETERMINATION PROCESS 


2.1 Introduction 

The attitude (state) of the balloon system can be determined as a 
function of time if (a) a method for simulating the motion of the system 
is available and (b) the initial state is known. The system motion can 
be simulated once the system model is determined. The initial state can 
then be obtained by fitting the system motion (as measured by sensors) to 
the corresponding output predicted by the mathematical model. In the case 
of the LACATE experiment the sensors consisted of three orthogonally 
oriented rate gyros and a magnetometer all mounted on the research plat- 
form. The initial state was obtained by fitting the angular velocity 
components measured with the gyros to the corresponding values obtained 

frxM the solution of the math model. 

A block diagram illustrating the attitude determination process 
employed for the LACATE experiment is shown in figure 2.1-1. The process 
consists of three essential parts; i.e., a process for simulating the balloon 
system (block 1), an instrxunentation system (block 2) for measuring the 
output, and a parameter estimation process (block 3) for systematically 
and efficiently solving the initial state. A more detailed discussion of 
each of these parts is presented below. 

2.2 System Simulation Process 

The main steps in the system simulation process are shown in the block 
diagram of figure 2.2-1. They consist of (a) development of a system 
mcdel to predict state (block 1) , (b)' solution of the model (block 2) and 
(c) developing a math model for computing the angular velocity conponents 
of the research platfom. Each of these steps will be discussed in greater 


detail below. 












vSystem math nKxjel A nystem math model must be obtained which enables 
one to predict the motion of tlie system at float alti.tude since all forms 
of output (c.g. angular velocity, tension in the cables, etc.) can be 
determined once this is known. This motion is very complex and involves 
various types of oscillation including bounce, pendulation and spin. 

Moreover, the complexity of the motion is increased with increasing number 
of subsystems. 

The math model of the balloon system is complicated primarily by the 

two important factors. They are as follows: 

(a) the balloon itself is a distributed parameter system which has 
motion in an infinite fluid media. Hence, it is necessary to first 
idealize it as an equivalent rigid body in order to develop a lumped parameter 
model for the entire system. 

(b) The balloon system is subjected to nondeterministic wind gusts 
which result in forces acting externally on the system. At the present 
time very little is known about the nature of these gusts. 

The exact dynamic model for the balloon itself consists of the 
equations of motion for the solid (i.e. balloon fabric) and the fluid 
dynamic equations. These equations are coupled through the boundary condi- 
tions which must be satisfied at the interface of the solid and fluid 
media. The resulting model is extremely complex and consists of a system 
of coupled partial differential equations. The system model is sijiplified 
by treating the balloon as a lumped parameter (rigid body) element. This 
is accomplished by developing approximate expressions for the afirodynamic 
forces and torques whidi result due to the interaction between the balloon 
and fluid media. These forces and torques are then treated as external 
reactions on the solid system. 

A linear systems model which includes the effect of the aerodynamic 
reactions has been developed (3) by neglecting the effect of second order 


terms. An equivalent form of this model which involves only the pendulation 
angles in two oi'thogonally oriented planes is given as follows 


0 + A0 ** 0 f 
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The elements are defined in terms of the system parameters in reference 
(3). The pendulation angles 0^ and and the .spin angle ^ -are illustrated 
in figure 2.2-2 and figure 2.2-3. 
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Solution of System Model . The solution to equation 2.2-1 can be obtained 
by assuming the following fozro for 0 ; i.o. 

0 = X sin nt. 2.2-7 

Substitution in to equation 2.2-1 yields the following eigenvalue problem, 

i.e. 

AX = fi^X, where 2.2-8 

^li 

X^(i = 1, 2, 3) = X 2 ^ are the eigenvectors and 

^3i 

JJ^Ci = 1,2,3) are the eigenvalues. 

Equation 2.2-8 can be solved numerically on the computer. A listing of the 
fortran program enployed to solve this equation for the particular problem 
under study is given in appendix A. 

The solution to equation 2 . 2-2 is obtained in the same way and yields 
the identical eigenvalue problem. The solution for equation 2.2-3 is 
easily obtained by integrating the equation twice with respect to time (t) . 
The final closed form solution for 0 , $ and ij >3 is given as 

0 (t) = ^l"^^ 

+ 5^(03 sin n 2 ‘t cos 02 "^) 

+ 5<2(Cg sin flgt + c 


Ig cos JJgt) , 


2.2-9 


^i(t) = °8 ^1^^ 


+ ^^(Og sin ^2^ <^X0 °°® ^2^^ 

+ sin fijt + Cj ^2 °o® ^3^5 ' 2.2-10 

<|>2(t) = at + (j>2(to) , where 2.2-11 

Cj (j « 1 , 2 ,..., 12) * unknown constants which are determined 

from the initial states; i.e. 0(to), 

• • 

51 (to), 0{to) and !jJ(to), 
a ® constant rate of spin, emd 
(j)^(to) = initial spin displacement. 

Computation Balloon Angular Velocities . The relationship between the 
platform motion (0^, S^, <l>^ and and the platform angular 

velocity components is obtained through the application of the Euler emgle 
transformations to the system platform shown in figure 2.2-4. The trems for- 
mation equations are given as follows; 

. • 

‘**2 " "^3 ’^3 *^3 *^3 "^3 ' 

U3 « 63 sin '1*3 + ^3 » where 2.2-12 

w^(i a 1,2,3) a angular velocity components of the platform 
along direction, 

63 a pendulation angle in the 63 ©3 plane, 

1II3 a pendulation angle in the e^ 63 plane, and 
a spin angle about the e, axis. 




The main object of the parameter estimation method is to determine 
the initial system state (Qq' ‘*’ 3 ° ^ rates ( w ) 

obtained from the rate gyros fit (over some time interval 0 ^ t ^ T) , in 
an optiml sense, those rates ( w ) predicted from the system model. With 
this initial state determined, the instantaneous system state is obtained 
sijiply as the output ( 0(t), if(t) and •t'jCt)) of the system model. Hence, 
the problem is basically one of the parameter determination in which the 
initial state parameters ( 6 © > ‘^' 3 ° ^ role of the unknown 

parameters. For the purpose of this work, the platform rates will be fit 
in a least square sense; i.e., a performance function ( $) will be formed 
and the initial state determined such that this function is minimized. 

The process will be repeated (i.e., 6 ^ , and (ji^o will he updated) 
every T seconds. Initially T will be set equal to 30 seconds (time required 
to acquire a complete radiance profile) although a study (to be conducted 
in the future) will be made to determine the minimum T required such that 
the necessary precision and accuracy are satisfied. 

In this research, the function $ is formed as follows; 

^ = X Z - S'. ) , where 2.3-1 

j=l i=l ^ 

N = number of data points taken in 0 ^ t _< T, 

s angular velocity conputed from system model, and 
= angular velocity given by the rate gyros and to be conpared with 

The function ♦ is clearly dependent on the initial state. This initial 
state is obtained from the condition that $ take on a minimum; i.e., by sol- 
ving the following optimization problem 

min. $ = $(w). 2.3-2 

The angular velocity conponents (u^) are obtained from the transformation 



equations 2.2-12 and since 0 , ^ ,(j)^ are functions of the initial state 
( 0 ^, and )» equation 2.3-2 can be written as 

min $ = $(0Q, ip^, (p^o) 2.3-3 

rw 

Since the values of $ are obtained numerically (o is given as a discrete 
data point) it will be necessary to employ some direct search technique to 
solve the above optimization problem. In general, the algorithm for any 
direct search techniques is given as follows 

- - k - k 

Xo = + ^^o ’ (k = 1,2,3,...), where 2.3-4 

is the vector of ddf .values , 

■■ ]c 

6X_ is a vector of increments, and 
o 

- k+1 

X^ is the vector of improved values. 

The vector SR^ is found such that 4 >(X^ ) < ) . The value of k is 

incremented until 5?^ converges; i.e., until the nom of satisfies sane 
error criteria. Figure 2.3-1 illustrates the application of equation 2.3-4 
to minimize $ . 

There are many direct search techniques for systematically determining 
6 After comparing several of these, the direct search method of Hooke 
and Jeeves was chosen for the following reasons, i.e., (a) the special 

feature in accelerating of distance, so called pattern search, and (b) 
this method is already available in subroutine form (5) and was errployed 
with the process outlined in figure 2.3-1. The discussion of Hooke and 
Jeeves direct search algorithm is presented in section 2.4. 

2,4 Hooke and Jeeves - Direct Search Method . 

The Ilooke and Jeeves [2] method is a sinple and powerful uni variant 
method for finding the minimum of a function. A modification of the basic 
univariant numerical search method, it involves trial explorations and then 














ever-expanding steps, called pattern nt>ves, in the direction indicated by 
the explorations. Tlie method is designed to follow a ridge, so that it does 
not suffer from the main disadvantages of the basic univariant approach 
which can not effectively cope with ridges or sharp valleys. 

The Hooke and Jeeves algorithm consists of two major phases, an 
"exploratory search" around the current base point and a "pattern search 
in the direction selected for minimization. Figure 2.4-1 is a simplified 
informational flow diagram for the algorithm as inplemented by \'Iook [43. 

The steps (blocks) in this figure are as follow: 

Block 1 The initial estimates for all decision variables (X^. ) as well 
as initial incremental changes or step sizes (AX^) in the decision variables 
are provided. 

Block 2 TTie objective function, ^>(X) is evaluated at the base point 
(X) which is the vector of initial guesses (X^) of the decision variables. 

Block 3 An exploratory search (type 1) is performed next, i.e. 
each decision variable (X^) is changed in rotation, one at a time, by the 
incremental amount (AX^) and the objective function is evaluated at the 
new point. If this incremental fails to improve the objective function 
then X. is changed by (-^X^) and the value of 0<) again evaluated as 
before. If the objective function is still not improved then X^ is left 
unchanged, and the same procedure is employed again with This 

process is repeated until all the decision variables (X^^, X 2 ,...,Xj^) haven 
been so changed. Figure 2.4-2 illustrates the steps ir. an exploratory search 
for a two dimensional problem. For each change in the decision variable, 
the value of the objective function # (R) is ccmpared with its value at the 
previous point. If, upon completion of the exploratory search, none of the 
changes yields an improved ^ (5?) (i.e. Xj^,X 2 ,. . • »Xj^) remain unchanged , then 
the stages in block 7 are performed next; otheivise block 4 is implemented. 

Block 4 /if ter completing the type I exploratory search, the new base 
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Figure 2 . 4-2 ' ExplQratoar/ Search Strate< 



point is set equal to the final base point obtained frxMii block 3. 

Block 5 After completing the type I exploratory search and obtaining 
a new base point, a "pattern search" is made. The new value of the decision 
variables define a vector, S(see figure 2.4-2), that represent a successful 
direction for minimization. A series of pattern searches is now made 
along this vector, usually in increments of 2|s| until $(R) no longer 
decreases. The magnitude of the step sizes for the pattern search (i.e. 

S'^ in figure 2.4-3) is roughly proportional to the number of successes 
previously encountered in each coordinate direction during-.the e>q)loratory 
searches for the previous cycle. The success or failure of a pattern 
move is not established until after a type II exploratory search (block 6) 
has been completed. If 4'(X) does not decrease after the type II explora- 
tory search, then the pattern search has failed and a new type I exploratory 
search (block - 3) is made in order to define a new successful direction. 

The base point for a pattern search in the case of a two dimensional problem 
is illustrated in figure 2.4-3. 

Block 6 The type II exploratory search (figure 2.4-3) is made after 
a tenporary exploration point is obtained fran a pattern search. The 
chief difference between the type I and typ>e II exploratory searches is the 
magnitude of the step sizes (iX^). In the case of the type II exploratory 
search, AX^ is taken as scm»e multiple of theAX^ (i.e. C/(aX£), see figure 
2.4-3) used in the type I exploratory search. This is done in order to 
accelerate the search. 

Block 7 If the type I exploratory search fails to give a new success- 
ful direction, then the current |ax^| is conpared to some preset allowable 
tolerance (error input in block 1). If |ax^| is larger than the allowable 
error, then block 8 is nrplemented. Failure to irirrove $ (X) for |ax^I 
smaller than the allowable error indicates that a local optimum has been 


r : Temporary base points 

(D , : Current base points 

® * Pattern search 



Figure 2*4-3 General Search Strategy 



reached and the search is terminated# 


Block 8 If I AX. I is still larger than the prespecified error, then 

' i ' 

I Ax^ I is reduced gradually until the type I exploratory seaix:h can be em* 
ployed to define a new successful direction. 

In order to terminate the search, two additional basic tests must 
satisfied, 'fhese are described as follows: 

(i) ,‘ifter each exploratory and pattern search, the increment in the 
objective function |a 4>| is compared with a prescribed test value. If this 
increment is less than the test value, then the exploratory of pattern 
search is said to have failed. In this case block 3 or 7 is ijipleirented. 

(ii) If |a 4>| is greater than the prescribed test value, then a test 
is made to determine if the objective has increased (a failures) or decreased 
(a successful search). This second test ensures that the values of the 
objective function is always being improved. 

The fortran coding for the Hooke and Jeeves method has been provided 
(with some minor revision) by M.I.T. Joint Conputer Facility. This is given 
in appendix B . The program is available in subroutine foirn with various 
parameters in the calling sequence. These parameters include 

(i) the number of decision variables (n) , 

(ii) the initial guesses of the variables 

(iii) the ..'equired accuracy (err), and 

(iv) the allowable number of allowable iterations. 

2.5 Test Problem 

In order to verify the accuracy and precision of proposed parameter 
estnjiHtion process (refer section 2.3) it is necessary to employ the process 
to a test problem in which the initial state is already known. Preferably, 
the test problem model should be identical to the one employed in this 
study; i.e. equatjons 2.2-1 to 2.2-3. With the initial state known, the 


resulting values of'o , ij; , (ji , 0 , ij/ and L which determines the 

•i *4 4 hS 4 

State of the platfom (as a function of time) can be obtained from equations 
2.2-9 to 2.2-11. With these, the values of uj and uij can be computod 
from equation 2.2-12. 

The testing process will consist of the following steps, 

(a) Assign fj.xed values to the initial state parameters 1§(0), JjKO), 

• • 

0 ( 0 ) and $(0), or equivalently to the unknown constants ’°12 ^ 

equations 2.2-9 and 2.2-10. 

(b) Confute the orientation (state) of the platform 0 ^, 6^, 1 ( 1 ^, 'l ’3 

and as a fvinction of time by employing equations 2.2-9 to 2.2-11. 

(c) Determine the values of ^2 and a function of time by 

substituting the results from part (b) into equation 2.2-12. 

(d) Etiploy the results (sanpled at various times) from part (c) as 
input to the parameter estimation process and utilize this process to 
recover the initial state values (or unknown constants Cj^, C 2 ,..., 

Tlie accuracy and precision of the process can be determined by com- 
paring the results of part (d) to the corresponding assumed values of part 
(a) . The testing prxxedura can be repeated for various sets of input 
values. A flow chart illustrating the test process is given in figure 2.5-1. 










CHAPTER III 


1 

?c, 


RESULTS AND CONCLUSIONS 
3.1 Data for LACATE Mission 

Ballcx>n and System Data . Figure 1.1-1 illustrates the actual I^acate 
balloon system and figure 1.1-2 illustrates the corresponding idealized 
system used in this study. The values for the various lengths and masses 
of the idealized system are given in table 3.1-1. 

The actual design profile of the balloon at float altitude is 
illustrated in figure 3.1-1 and the corresponding dimensions are given 
in table 3.1-2. Values for (a) mass center, (b) center of volume, and 
(c) mcmient of inertia (at the mass center) of the balloon were confuted and 
these values are presented in table 3.1-3 along with other balloon pro- 
perties . 

Standard Atinosphere Data . Figure 3.1-2 illustrates the U.S. Standard 
Atmosphere and the corresponding properties are presented in table 3.1-4. 
Table 3.1-5 presents the actual properties of the atmosphere at the float 
altitude (47 kilometers). 

Math Model Data . The expressions for the element of A matrix of 
equations 2.2-1 and 2.2-2 for the system model are given in refe.rence (3). 
The numerical values for these elements were conputed based on the data 
given above and these values are presented in table 3.1-6. 

Gyro Data . The platfom coordinate axis which were enployed for 
referencing the angular velocity components (as measured by the instrumen- 
tation packcige) do not coincide with the coordinate axis used in the Euler 
angle transfonnation equation (2.2-12). The relationship between these two 
coordinate systems is shown in figure 3.1-3. The resulting transformation 
equations which relate the angular velocity components are given as 
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«1 S 

> 

3.1-1 

«2 = '^2 

, and 

3.1-2 

= -■S?3 

, where 

3.1-3 


i»£ is the angular velocity component along the 
axis, and 

is the corresponding angular velocity conponent 
obtained from the gyro. 

The numerical values of the angular velocity components obtained frcmi 
the gyro system are presented in appendix C. Typiccil plots of the angular 
velocity conponents obtained from the gyros are illustrated in figures 
3.1-4 to 3.1-6. 

The azimuth angle obtained from the magnetometer was measured 
clockwise from magnetic north to the negative platform axis. However, 
the angle ^ in equation 2.2-11 is measured counter-clockwise fron north to 
the same negative X^^ platform axis. This is illustrated in figure 3.1-7. 

The transformation equation relating these two angles is given as 

(|> s 360° - 4 , where 3.1-4 

m 

^ is the spin displacement in equation 2.2-11, and 

(f>j^ is the spin displacement measured by the magnetoneter. 

A typical plot of the azimuth angle <^^is shown in figure 3.1-8. 

Eigenvalue Problem . The solution to the eigenvalue problem (equation 2.2-8) 
was obtained by enploying the ccnputer program given in appendix A in con- 
junction wi1di the coefficients pres&iited in table 3.1-6. The solution for 
the eigenvalues and corresponding eigenvectors is presented in table 3.1-7. 
The magnitude of 0^ is equal to the natural frequency of the system. 


TABLE 3.1-1 


Idealized LACATE System Properties 


m3 


(maximum radius of the balloon) = 249.989 ft. 

(distance fron center of mass of balloon shell to mass m^^) = 222.091 ft, 
(distance from mass m^^ to m 2 ) = 75 ft. 

(distance from mass m 2 to m 3 ) = 15 ft. 

(distance fnan center of mass (C.G.) to center of volume (C.V.)) = 23.712 ft. 
(mass of balloon shell) = 2850 Ib^^^ 

(lumped mass) = 135 Ib^^ 

(lumped mass) = 135 Ib^^ 

(lumped mass) = 375 Ib^^^ 












TABLE 3.1-3 


BALLOON PROPERTIES 

(maximum radius) = 249.989 ft. 

\ (height, corresponded to r^) = 209.526 ft. 

Vjj (inflated volume) = 45,378,282 cu. ft. 

H (inflated height) = 350.1 ft. 

^o weight of balloon shell including top cap weight) = 2880 lbs. 

Gore length = 674.83 ft. 

Surface area = 635,711 sq. ft. 

(distance from bottom apex to the mass center of balloon shell) = 222.091 
^2 (distance from nadir to center of volume) = 198.379 ft. 
t (thickness of balloon shell (strato film®) = 0.0006 inch. 

(moment of inertia at the mass center (C.G.) = 2.622374 X 10^ slug-ft^ 
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Table 3. 1-4 

General Properties of the U.S. Staindard Atmosphere (6) 


ALTI- 

TUDE, 

(m) 

TEMPER- 

ATURE 

(C) 

TYPE OF 
ATMOSPHERE 

LAPSE 

RATE 

(C/km) 

‘‘ 1 
(m/scc^) 

A 

PRESSURE 

DENSITY 

p(kg/m*) 

e 

15.0 





L013 X 

10^ 

1.225 


- 


Potylropic 

-6.5 

9.790 

1.235 





11,000 

-56.5 





226i K 

10^ 

3.639 

X lO"*- 



Isothermal 

0.0 

9.759 






20,000 

-56.5 





5.475 X 

10^ 

8.804 

X 10"* 



Poly tropic 

+ 1.0 

9,727 

0.972 





32,000. 

-44.5 





8.680 X 

10^ - 

1.323 

X 10"* 



Polytropic 

+2.8 

9.685 

0.924 





47,000 

-2.5* 





1.109 X 

10^ 

1.427 

X 10"* 



Isothermal * • 

0.0 

9.654 


* 


• 
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52.000 

-2.5 
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lO’ 

7.594 

X 10"* 



Polytropic 

-2.0 

9.633 

L0£I 
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-20.5 





I.82I X 

10* 

2.511 

X 10"^ 



Polytropic 

-4.0 

9.592 

U36 





79,000 . 

-92.5 





1.035 


2.001 

X 10"* 


- 

Isothermal 

0.0 

9.549 






88,743 

-92.5 





1.644 X 

I0~‘ 

3.170 

X 10"‘ 


C = temperature in degrees Centigrade 
km = Killometers 
n = polytropic exponent 
g = local acceleration of gravity 


TABLE 3.1-5 


Prdp^ies of 'Atifbsphei^e at float Altitude (47 'Km) 

pQ^ (density of air) = 2.78 X 10“^ 

(density of helium) s 3.834 X lO""^ 

Pq (pressvtre) = 2.254 ^2 

Cjj (viscous drag coefficient) = .5 

g (gravity) = 31.7 ^^2 

(viscosity of air) = 3.57 X lO""^ ^ ^ 
y„ (viscosity of helium) = 3.74 X 10“^ 

ii XX* S6C • 


I 


TABL£ 3.1-6 

Pr^Amelents of A Matrix [Equation 2.2-1 8 2,2^ 


^11 ’ 
^21 ' 
^31 
^12 " 
®22 " 
^32 “ 


2.10463 X 10 ‘ 
-6.08078 X 10'^ 
-5.07147 X 10"^ 
-3.85241 X 10"^ 
1.74685 
-7.98370 


a23 = -1.17407 
a33= 7.98370 



= coordinate system enployed for purposes of math model. 
X^Q s gyro coordinate system. 


Figure 3.1-3 


Coordinate System 













imuth Anqle (<^ ) From Magnetometer 
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3,2 Ev aluation of Parameter Deterinination Precess 

The process for evaluating the attitude detemination process by 
eirploying a test problem was discussed in section 2,5. The form of the 
jiBth nodel for the test problem was identical to the system model given in 
equation 2,2-1 through 2.2-3. Several different sets of unknown constants 
(C.) were assumed and the resulting angular velocity components (oi. ) were 
computed by employing the solution equations 2.2-9 to 2.2-11 and the trans- 
formation equation 2 . 2-12 . These values of oJ were input to parameter 
determination process which was then employed to recover the unknown con- 
stants . 

Table 3.2-1 gives the results for one test case. These results indi- 
cate that the process is capable of recovering the unknown constants 0 ^^^ 
with good precision. The results presented in table 3.2-1 are typical of 
those obtained for other test cases. Hence one can conclude that the process 
is suitable for evaluating the attitude of the LACATE S 3 ^tem. 

During the next stage of the research, the process will be employed 
(in conjunction with the math model given in eqs. 2.2-1 - 2.2-3) to deter- 
mine the angular velocity components of the platfonn of the LACATE balloon 
system. These values will then be oompared to the actual data in order to 
evaluate the system model. Results will be presented in the next report. 
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